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We performed first-principles calculations to investigate the possible magnetism induced by the 
different concentrations of non-magnetic impurities and vacancies in BN sheet. The atoms of Be, B, 
C, N, O, Al and Si are used to replace either B or N in the systems as impurities. We discussed the 
changes in density of states as well as the extent of the spatial distributions of the defect states, the 
possible formation of magnetic moments, the magnitude of the magnetization energies and finally 
the exchange energies due to the presence of these defects, ft is shown that the magnetization 
energies tend to increase as the concentrations of the defects decreases in most of the defect systems 
which implies a definite preference of finite magnetic moments. The calculated exchange energies 
are in general tiny but not completely insignificant for two of the studied defect systems, i.e. one 
with O impurities for N and the other with B vacancies. 

PACS numbers: 75.75.+a, 73.22.-f, 67.57.Pq, 61.72.Ji 



I. INTRODUCTION 



The magnetism involving only s- and p-electron ele- 
ments continues to attract much attention due to the 
potential of extensive applications as well as the urges to 
understand its physical origins. Recently, some experi- 
mental groups have discovered either weak or strong fer- 
romagnetism in fullerenesp], H and graphite [3[ systems. 
A few theoretical studies attempting to find magnetism 
in some potential non-magnetic systems have also been 
carried out previously However, the origin of ferro- 
magnetism in those systems is still under debate on both 
experimental and theoretical sides [3, The mecha- 
nism for forming magnetic ordering in solids, such as fer- 
romagnetism, antiferromagnetism, etc., can be referred 
to the direct or indirect exchange interactions (super- 
exchange, double exchange, etc.) among magnetic mo- 
ments. Therefore, a strong enough exchange interaction 
should be the crucial criterion to determine this possi- 
ble new class of magnetic materials involving s- and p- 
electron atoms. In this work, the boron nitride(BN) sheet 
is used as the host system and the possible formation 
of magnetic moments by non-magnetic defects including 
substitutions of impurity atoms and creations of vacan- 
cies is studied. The exchange energies of those systems 
possessing finite magnetic moments are further examined 
and we shall demonstrate that there exist mostly very 
weak interactions among these defect-induced magnetic 
moments except two cases. 

BN can form three different bulk structures which 
are hexagonal BN(h-BN), cubic BN(c-BN) and wurtzite 
BN(w-BN). Of these three structures, h-BN is the room 
temperature phase. Similar to graphite, h-BN is quasi- 
2D with weak interaction between layers. Nevertheless, 
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different from the delocalized n electrons in graphite, the 
stronger electronegativity of N than that of B causes it 
electrons to distribute more around N . This strong direc- 
tional effect of bonding confine the motion of 7r-electrons 
and thus results in a gap in h-BN. BN is the lightest III- 
V compound of those that are isoelectronic with III-V 
semiconductors such as GaAs, but with wider band gap 
i.e. Eg(BN) = 4.0 ~ 5.8 eV at room temperature^, 
compared with 1.42 eV of GaAs. 

A 2D-structural material is usually a good basis for 
studying the physical properties of the most interesting 
nanostructures, such as nanotube, nano-ribbon etc.. The 
tubular structure of BN had been synthesized [7| experi- 
mentally. Moreover, it has been shownQ that, in differ- 
ent ways of rolling up the BN nanotube, the electronic 
properties are all semiconducting which is the same as the 
extended BN sheet. This is very different from that of 
graphite and carbon nanotube which can be either metal- 
lic or semiconducting, depending on the way of rolling up. 
Therefore, one would expect that the effects of defects in 
BN sheet can catch many important physical properties 
as those would exist in the BN nanotubes. 

In semiconducting and insulating systems, magnetic 
moments can be induced by defects [3. Hoi fl~TT| . Here, two 
different types of defects are used to study the possible 
formation of magnetic moments in BN sheet. One is 
that the atoms of Be, B, C, N, O, Al and Si are used 
to replace either B or N as impurities. The other is to 
create vacancies by removing either B or N atoms in BN 
sheet. 

Our plans for this article are as the followings: Sec. 
II summarizes the computational methods in this work. 
The results of non-spin-polarized calculations are pre- 
sented in Sec. III. In this section, the Subsec. A will 
present the electronic properties of BN sheet by the an- 
alyzed orbital-projected partial density of states (DOS). 
The Subsec. B includes the effects on DOS due to de- 
fects, i.e. the location of the defect states in DOS, as 
well as the extent of the spatial distribution of these de- 
fect states. All the results of spin-polarized calculations 




FIG. 1: (Color online) The 4x4 unit cell for the hexagonal 
BN sheet. D represents defects created in the system, i.e. 
impurities or vacancies. 

are presented in Sec. IV. In the Subsec. A, we will inves- 
tigate the possible formation of finite magnetic moments 
induced by defects while the results of exchange energies 
and the corresponding Curie temperatures will be pre- 
sented in the Subsec. B. Finally, the conclusions are in 
Sec. V. 



II. COMPUTATION METHOD 

Our calculations are based on density functional 
theory (DFT)[l2j with generalized gradient approxima- 
tion(GGA) of Perdew and Wang[13| for the exchange- 
correlation energy functional. The PAW method[l4| im- 
plemented by Kresse and Joubert[15] is used to describe 
the core- valence electrons interactions. There are 2, 3, 
4, 5, 6, 3 and 4 valence electrons included here for Be, 
B, C, N, O, Al and Si respectively. The one-electron 
Kohn-Sham wavefunctions are expanded by the plane- 
wave basis with kinetic energy cut-off (E cut hereafter) of 
400 eV. The sampled k-points in the Brillouin zone(BZ) 
are generated by the Monkhorst-Pack technique [l6[ with 
gamma centered grids for hexagonal lattice. All calcu- 
lations are carried out by Vienna ab initio simulation 
package(VASP)|l3. 

The calculated lattice constants of h-BN are a — 2.5A 
and c = 6.36A with 13 x 13 x 5 k-point grids and 
compared with the experimental results (a = 2. 5 A and 
c = 6.66A)pj| are reasonably well[l9j]. The supercells 
composing of 4 x 4 primitive unit cells of BN sheet are 
used to simulate systems with defects(See FIG. [1}. The 
vacuum distance between BN sheets (taken as along the z 
direction) was chosen to be around 13A which was tested 
with 3x3x1 k-point grids to be large enough to avoid 
interactions between the sheets. 

Of all the relaxed configurations in this work, 
the atomic forces calculated by Hellmann-Feynman 
theorem[2l| were smaller than 0.02 eV/A. In the spin- 
polarized calculations, we performed the above atomic re- 
laxations again to ensure the relaxed configurations. The 
DOS and magnetic moments of the final relaxed struc- 
tures were then evaluated. The larger supercells of 8 x 4 
and 8x8 were used to study the possible finite magnetic 
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FIG. 2: (Color online) The orbital-projected partial density 
of states for N and B atoms in hexagonal BN sheet. 



moments induced by defects, as well as the variation of 
the magnetization energy (denoted as Em hereafter) and 
the exchange energy J with respect to the distances be- 
tween defects, i.e. defect concentrations. 



III. NON-SPIN-POLARIZED CALCULATIONS: 
DENSITY OF STATES 

A. BN sheet 

The orbital-projected partial DOS of B and N in BN 
sheet are presented in FIG.[U The colors of blue, yellow, 
green and red in the figures represent the s-, p x -, p y - and 
p z -orbital-projected partial DOS respectively. The spiky 
feature of the DOS is due to the dimensional reduction 
from 3D bulk h-BN to 2D planar structure of BN sheet. 
Two valance bands (shorted as VB hereafter) and the 
lowest 3eV range of the conduction band (shorted as CB 
hereafter) are included in the figure for discussion. The 
DOS's in the same energy range for BN sheets with de- 
fects will be used in the following discussion for compar- 
ison and the bottom of the CB is taken as the zero of the 
energy for convenience in discussing the splitting of the 
electronic bands from the two VBs as well as the CB due 
to the presence of defects. The two VBs in BN sheet arc 
well-separated as those in the bulk h-BN[l8j]. The lower- 
energy VB (VB1 hereafter) contains 2 electrons: one is 
from the s-electron of N and the other is from the hy- 
bridized s- and p-electrons (sp 2 -electrons hereafter) of B. 
The higher-energy VB (VB2 hereafter) contains 6 elec- 
trons, i.e. three sp 2 -electrons and one Pz-electron from N 
and two sp 2 -electrons from B. From FIG Owe see that 
the nature of the two VBs for B is mainly hybridized 
suborbitals while that for ./V is s-like orbitals in VB1 and 
hybridized sp 2 -like orbitals in VB2. Near the top of VB2 
and the bottom of CB, the DOS's are mostly p 2 -orbital 
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FIG. 3: (Color online) The DOS's for the BN sheet with 
defects. The upper pannel shows defects on B while the lower 
one shows defects on N. (a)-(f) are the results of the defect 
systems with impurities Be, N, B, C, O, Al and Si atoms 
respectively. The defect systems with B and N vacancies are 
shown in (g). All dashed lines cut at Bf- 



in nature. However, the dominant contribution switches 
from the p z -orbital on N atoms near the top of VB2 to 
the ^-orbital on B atoms near the bottom of CB. That 
is the electronic excitation in BN sheet would involves 
displacement of electron distribution from the p 2 -orbital 
around N atoms to the p z -orbital around B atoms. 



B. Effect of defects 

In the present section, we discuss the effect of defects 
on the DOS of BN sheet when the spin-polarization is 
not yet switched on. We shall demonstrate that the ef- 



FIG. 4: (Color online) Schematic diagrams of the effect of 
defects on DOS. The dashed lines show the possible circum- 
stances in forming metal-like DOS. 



feet of defects on the DOS can be generally categorized 
into two groups according to the chemical properties of 
the defects. The 4x4 supercell was employed in all the 
results presented in this section which corresponds to a 
doping of 1/32 (~ 3%) of the defects in the system and a 
distance of 10.03A between the neighboring defects. The 
possible interaction between defects will be discussed in 
the following section when the exchange energy is inves- 
tigated. However, the results presented in the present 
section are applicable to the systems with low-density 
dopings of defects as we shall demonstrate in the follow- 
ing section that the interaction between defects is already 
quite small in the present configurations. The effect of 
C impurity on the DOS of BN sheet will be discussed in 
details first and then followed by the results for the series 
of defects we studied. For simplicity, we denote, e.g. Cb, 
as the defect system which consists of an impurity of C 
atom substituting one of B atoms in the original super- 
cell of BN sheet and Vb(Vn) as the system consisting of 
a vacancy replacing one of the original B(N) atoms. 

One of the effects of the impurities is to introduce elec- 
tronic states into the energy regions which originally have 
no DOS in the BN sheet. We present the total DOS of 
all the defect systems considered in this study in FIG [3] 
In the case of Cb (FIGO^c)), the defect states split from 
the original three bands and move to the lower energy 
region. This is due to the stronger bonding ability of 
C ions for electrons compared to that of B ions. On the 
contrary, in the case of Cjq , the bonding ability of C ion is 
weaker than that of N ion for electrons such that the de- 
fect bands move to the higher-energy positions compared 
to the original three bands. Therefore, the locations of 
defect bands is determined by the bonding ability of the 
defects compared to the substituted host atoms. The ef- 
fects of defects on the DOS's can be summarized into two 
types as outlined schematically in FIG(U For Nb, Cb, 
Ob, On and Sis, the bonding abilities of the impuri- 
ties are stronger than the substituted host atoms. We 
name these as type I. For Bes, Bejy, -Bat, Cn, AIn and 
Sijsi, which are named as type II, the three defect bands 
all move to the region of higher energy than the original 
three bands. 

Of these studied defect systems, some consist of odd 



4 



number of electrons in the supercell, i.e. Bes, BeN, Cb, 
Cn, Ob, On, Sis, Sipf, Vg and Vat, such that there is 
an unpaired electron occupying the defect state located 
between the VB2 and CB for BeN, Cb, On, Ob, Six 
and Vn, the top of VB2 for Bes and Vb as well as the 
bottom of CB for On and Sib- We denote DB as the 
state that the unpaired electron occupies hereafter. We 
shall show in the next section that magnetic moments are 
likely to form for those cases with partially occupied DBs. 
In case the partially occupied defect bands which form 
definite magnetic moments are also extended in nature, 
the long-range magnetic order can develop. Hence, it 
is expected that for Bes, On, Sis and Vb whose defect 
states locating at the edge of the original extended bands, 
the defect bands are more likely to be extended in nature 
and then more likely to lead to larger exchange energies 
as will be presented in the next section. 

In order to determine the properties of the partially- 
filled DBs, the extent of the states in space is established 
from summing over the contributions of the orbital- 
projected partial DOS of the DBs for the neighboring 
atoms of the defects as presented in FIG[5] For Bes, 
BeN, Ob and Oat, we see that the electrons of the DBs 
distribute more on the nearest neighboring atoms of the 
impurity than the impurity itself. This is contrary to the 
Cb, Cn, Sib and Sin systems in which most of the elec- 
trons of the DBs distribute around the impurity atoms. 
These results demonstrate that the electrons of the defect 
states due to impurities do not have to distribute mostly 
on the impurities themselves. About the spatial extent 
of the defect states, FIGf5]shows that only three systems, 
i.e. Bes, On and Vn, reach beyond 5 A. However, the 
defect states for all the defect systems considered in the 
present study are confined in a distance of 6A around 
the defects no matter whether the partially-filled defect 
states are at the edge of the extended bands or not. Fi- 
nally, we should emphasize that FIG [5] also reveals two 
categories of DBs: one is formed by the sp 2 -hybridized 
orbitals, e.g. s, p x and p y , and the other is p z in na- 
ture. These orbital-projected-DOS analyses imply the 
difference of the electronic distributions of the magnetic 
moments (detail discussions in the next section) formed 
by the unpaired electrons in DB, e.g. the moments are 
distributed on the plane of BN sheet for BeN, Ob and 
Vb while the moments are distributed perpendicularly to 
the plane of BN sheet for BeB, Cb, Cn, On, Sib, Sin 
and Vn systems. 
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FIG. 5: (Color online) The ratio of the orbital-projected par- 
tial DOS summing over the ith nearest neighbours of the de- 
fect of the system with partially-filled DB. There are l(the 
defect), 3, 6, 3, 6, 6, 3, 3 and 1 atoms contributing to each 
bars in sequence, a is the lattice constant 2.5A. 



A. The possible formation of magnetic moments 



IV. SPIN-POLARIZED CALCULATIONS: THE 
EFFECTS OF DEFECTS 

In this section, we present the results of spin-polarized 
calculations to investigate the possible formation of finite 
magnetic moment in the BN sheets containing defects as 
well as the possible long-range magnetic ordering in these 
systems. 



To identify the possible magnetic moments due to de- 
fects as well as the extent of the magnetic moments, the 
4x4 and 8x8 supercells consisting of only one defect 
were employed to simulate the BN sheet with defect con- 
centrations of 3.125% and 0.78% respectively, i.e. cor- 
responding to a distance of 10.03A and 20.05A between 
nearest-neighbour defects (denoted by do hereafter). In 
order to obtain the relevant physical quantities for the 
cases within the range between do — 1QA ~ 20 A, we 
create a rectangular supercell as shown in FIG1H1 By 
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FIG. 6: (Color online) The rectangular supercell with cell 
dimensions 20.05A and 17.37A. The distance between the 
two defect sites is 13.26A. The supercell contains 64 primitive 
cells of BN sheet. 
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FIG. 7: (Color online) The calculated magnetic moments with 
respect to do, the distance between the nearest-neighboring 
defects. 
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FIG. 8: (Color online) The magnetization energy Em v.s. di 



FIG. 9: (Color online) The relaxed structure for the defect 
system of Ob- The impurity, O atom, was relaxed and moved 
toward two of the three neighboring N atoms to open up one 
of the original N — O bonds. The spin-polarized energy for 
this configuration with unequal N — O bond lengths is lower 
than that with equal N — O bond lengths by lllmeV. 



providing cither one defect or both defects in the super- 
cell of FIGO one can simulate the cases of 17.37A and 
13.26A for dg respectively. The calculated magnetic mo- 
ments of the systems with respect to dg are presented 
in FIG 13 For Ng, Bn, Alg and AIm which have even 
number of electrons in the supercell and completely-filled 
DBs, there exist no magnetizations and they are therefore 
excluded from the following discussions. The calculated 
magnetic moments in Beg and Vg were found to vary 
as dg increases. However, their magnitudes converge to 
Ifig as dg increases to 17.37A. For Bex, Cg, CV, Og, 
On, Sig, Six and Vn, the calculated magnetic moments 
are always found to be \[ig and independent oidg. Note 
that for all the defect systems possessing finite magnetic 
moments, the defect bands are always partially occupied. 
The \[ig of magnetizations in these defect systems should 
be attributed to the unpaired electron occupying the de- 
fect states. 

To understand how stable these finite moments in BN 
sheet are, the magnetization energies, i.e. Em, obtained 
from the total-energy difference between the systems 
with and without spin-polarized configuration were cal- 
culated. The results are plotted in FIG[5] Firstly, the 
Ems of these systems tend to increase as dgs increase. 
This suggests a definite preference of finite magnetic mo- 
ments for these defect systems with distant non-magnetic 
defects. Secondly, the corresponding Em for these de- 
fect systems with distant defects can be determined by 
the saturated values in the figure, e.g. 360meV for Vg, 
300meV for Og, 220meV for Be^ etc., and then are 
found to be in the wide energy distribution (from 16meV 
of Beg to 360meV of Vg). We notice that the defect 
systems with metal- like electronic properties (Beg, On, 
Sig and Vg), i.e. with fermi level at the edge of either 
VB2 or CB, do not necessarily lead to large Em- For ex- 
ample, the Em of Beg is very small, i.e. 16meV which is 
even smaller than the room temperature thermal energy 
of 25meV. However, the Em of Vg is the largest, i.e. 
360meV, of all the systems we considered. We should 
mention that, when the defect distance is increased, the 
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metal-like DOS does not always hold. The bandwidth of 
the defect bands might become smaller due to localiza- 
tion such that they eventually separate from the original 
VB2 or CB bands. Of the studied systems with metal- 
like DOS at do = 10.03A, this happens in Sis and Vb 
whose defect bands are originally at the edge of CB and 
VB2 respectively. For Bes and On, the defect bands 
remain at the edge of VB2 and CB respectively up to the 
largest do = 20.05A we studied. 

Finally, for those systems with p z -\ike magnetic mo- 
ments, the bond lengths between defect and neighboring 
host atoms were found relaxed but preserving the three- 
fold symmetry. However, we found that, the relaxed 
bond lengths of BeN, Ob and Vb which have planar dis- 
tributions in the magnetic moments (the sp 2 -hybridized 
electron in DB) undergo structural distortions and break 
the three fold symmetry of the original BN systems after 
the spin-polarization calculation is switched on (John- 
Teller(JT) effect). FIG|5] shows the relaxed structure of 
Ob- The impurity, O atom, moved toward two of the 
neighboring N atoms to open up one of the three orig- 
inal N — O bonds. The total energy of this distorted 
structure with spin-polarization is lower than that of the 
structure with equal N — O bond length by lllmeV, simi- 
larly for Be n by 33meV and for Vb by 47meV. The effect 
of JT distortion are also responsible for the consequences 
of relatively larger magnetization energy of Bcn, Ob and 
Vb with planarly distributed moments compared to those 
with p z -like moments. 



B. The exchange energy 

From the above analyses, we have demonstrated that 
the non-magnetic defects can induce finite magnetic mo- 
ments in BN sheet. To identify the possible long-range 
magnetic orderin g in these systems, the Heisenberg-type 
of spin couplings [22j is employed to model the interac- 
tion (J, i.e. the exchange energy) between the nearest- 
neighbor magnetic moments due to defects. The ex- 
change energy J is determined from the total-energy dif- 
ference between two spin configurations, i.e. one with an- 
tiparallel and the other with parallel spin configuration 
whose energy are denoted as Eafm and Efm respec- 
tively hereafter. For simulating two defects distancing 
10.03A, we employed either the 8x4 supercell with two 
defect sites in each sub-supercells of 4 x 4 or the 8x8 
supercell with those in the two diagonal sub-supercells of 
4 x 4 to generate different spin configurations. In the anti- 
ferromagnetic configuration, the 8x4 supercell describes 
a configuration of six nearest neighbours (do = 10.03A) 
in which four are antiparallel- and two are parallel-spin 
to the centered defect while the 8x8 supercell describes a 
configuration of two antiparallel-spin nearest neighbours 
(d D = 10.032) and a d D = 17.37A for the four next- 
nearest neighbours. The rectangular supercell in FIG|5] 
with defects on both sites is also used to simulate the two 
spin configurations for the systems with defect distance of 
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FIG. 10: (Color online) The exchange energy 8J v.s. defect 
distances where J = ^{Eafm — Efm}- The energy range 
from — 5meV to 5meV is the estimated region where the mag- 
nitude of 8J is too small to be used to identify the sign of 8 J 
in our calculations. The dashed line indicates the value of ex- 
change energy respected to 50K of Curie temperature (T C MF ) 
under mean-field approximation. 

13.26 A. In the rectangular supercell, there are four near- 
est neighbours for a defect and all are antiparallel-spin 
neighbours in the antiferromagnetic configuration. We 
have used both 8x4 and 8x8 supercells for the Vb and 
On systems and the exchange energy obtained from the 
8x8 supercells were considered in the discussion as they 
correspond more to the low-doping conditions we would 
like to consider. The different values of the exchange 
energies obtained from using 8x4 and 8x8 supercells 
indicate that the simple pair-wise Ising interaction is not 
a good enough model to describe the energies of these 
systems of do = 10.03A We notice that both Vb and 
O n have DB at the edge of either CB or VB2. The 8x8 
supercells were also used on a few of the other defect 
systems, but all lead to similar results as those obtained 
from the 8x4 supercells. 

Our results are summarized in FIG I10I where the lines 
are to guide the eyes. From the previous magnetization 
energy calculations, our computational accuracy can rec- 
ognize meaningful energy difference of 5meV and larger 
values between configurations. For those systems of 
Bes, BeN, Cb, Ob and Sin, the calculated 8Js which 
are less than 5meV are too small to determine whether 
the systems are ferro- (8J > 0) or antiferro-magnetic 
(8 J < 0). For On, Vb and Sis, they are found to be fer- 
romagnetic, and for Cn and Vn, anti- ferromagnetic at 
do = 10.03A We should mention that the relaxed spin 
configurations of all the systems discussed here remained 
anti-ferromagnetic at the end of the calculations if with 
initial antiparallel-spin configurations. The exchange en- 
ergies for On, Vb and Sis, which possess metal-like 
DOS, are relatively larger at do = 10.03A, but not for 
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the case of Be b whose Em is tiny anyway. Although the 
magnetization energy of Vb can be as large as 200meV at 
do = 10.03A, the exchange energy is considerably much 
smaller (5.6meV), similarly for Sib whose Em ~ 90meV 
while J is smaller than 2meV. (Note FIG [10] shows the 
data of 8 times the exchange energy J.) This contrast 
is much reduced in the case oi On whose Em ~ 17meV 
and J ~ 9AmeV at d D = 10.03A. When the do is in- 
creased to 13.26A, the magnitude of exchange energies 
for all the above systems decrease. At the end, there 
are only two systems, i.e. On and Vb, which have large 
enough exchange energies to be identified numerically as 
the magnetically ordered systems. 

Within the framework of Heisenberg spin model, it is 
possible to estimate the Curie temperatures (T c ) of these 
systems. The mean-field result is given by 

k B T c MF = 2 - Jo, (1) 

where Jo is the on-site exchange parameter reflecting the 
exchange field created by all the neighbouring magnetic 
moments [H Hi] . Note that the T c determined by Eq.(l) 
is regarded as an overestimated values[25j]. Since the 
physics of T c involves the disorder nature of the system 
due to temperature, Jo should be related to the energy 
difference between the local disorder state and ferromag- 
netic state of the system[24]. However, the exchange en- 
ergies obtained here do not include disorder effect, i.e. 
Eafm is calculated from ordered anti-ferromagnetic con- 
figuration. Hence, the estimated T c (denoted as T^ F 
hereafter) by just simply substituting J in FIGQIJ] into 
the Jo in Eq.(l) is expected as an upper-bound of T C MF 
for a given system. For On, Vb and Sis with relatively 
large Js at d D = 10.031, the f c MF s are 72K, 42>K and 
13K respectively. However, they reduce to around 2QK 
for On and Vb, below 5K for Sis when do is increased 
to 13.26A In addition, the T™ F s of the defect systems 
with exchange energies smaller than 5meV are all below 



5K in different concentrations. Therefore, one can con- 
clude that the effect of the long-range magnetic ordering 
is very weak in all the defect systems considered here ex- 
cept for the defect systems of On and Vb- However, our 
calculations suggest that the systems of BeN, Cb, Cat, 
Ob, On, Sib, Sin, Vb and Vn all have finite magnetic 
moments and are therefore at least paramagnetic. 



V. CONCLUSIONS 

Different concentrations of non-magnetic impurities 
and vacancies in BN sheet were studied using first- 
principles methods to investigate the possible magnetism 
in these systems involving only s- and p-electron ele- 
ments. We firstly studied the effects of these defects on 
DOS and analyzed the characters as well as the spacial 
extent of these defect states. The magnetization energies, 
possible magnetic moments as well as the exchange en- 
ergies for these defect systems in different concentrations 
were evaluated. We demonstrated that all the defect sys- 
tems with partially-filled defect bands exhibited a definite 
preference for finite magnetic moments. The calculated 
exchange energies for low-density defect systems are all 
tiny except for the On and Vb whose exchange energies 
are not completely insignificant with an estimated T c of 
20K. 
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